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Abstract 

We consider the waveguiding by thin patterned slabs embedded in a homogeneous medium. In 
the longwave limit, the wave spectra of slabs are found to be well described by a single frequency- 
independent parameter, which we call the "guiding power". The guiding power can be evaluated 
in an effective medium approximation, similar to the Maxwell Garnett theory, but modified for 
the local field corrections specific to the two-dimensional geometry. The guiding power is different 
for the transverse magnetic (TM) and transverse electric (TE) polarizations. We show that the 
confinement factor of TM waves in a porous layer with high index ratio can exceed that for a 
homogeneous layer. Similarly enhanced confinement of TM waves is demonstrated for a layer of 
elongated cylinders or elliptic inclusion with a high axis length ratio. The effect originates from 
the suppression of local field effects and the increasing internal field in the inclusion. It may be 
useful in the design of far-infrared or THz quantum cascade lasers. 

PACS numbers: 42.70.Qs, 42.79.Gn, 41.20.Jb, 42.25.Lc, 78.66.Qn 
Keywords: waveguiding, modal control, photonic crystals 
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I. INTRODUCTION. GUIDING POWER 



Waveguiding of light in layered patterned structures, such as slabs of two-dimensionally 
periodic photonic crystals (2D PC), has attracted much interest in view of potential photonic 
applications The studied 2D patterns include periodic lattices of deep etched 

air pores or "conjugate" lattices of high permittivity cylinders. The patterned slabs (not 
necessarily periodic) can be employed in waveguides as either a core or a cladding. Numerous 
theoretical computations of the band spectra of 2D PC and PC slabs have been reported, 
based on expansions of the electromagnetic field in plane waves M 4| or cylindrical waves 
5, 6], as well as based on finite-difference time-domain methods {7, 8|. 

The low-frequency region of electromagnetic waves in the 2D PC is well understood. 
The waves have a linear spectrum that is very close to that obtained in the effective media 
appreciation flQ with an effective permittivity corresponding to the Maxwe.. Garnet, 
theory [9|,H1|. Thi s 

means that when the wavelength A exceeds the structure period a, the 
optical properties are primarily determined by the filling factor / of the inclusions (i.e. their 
total volume fraction) and do not rely on their long-range order or their shape variation. 
The disorder leads to a weak (for A ^> a) Rayleigh-like scattering. 

In this paper we investigate the waveguiding by PC slabs in a long wavelength frequency 
range, qd <C 1, where q is the wave vector, d is the thickness of the active or a core layer of 
the waveguide. We show that in this range the waveguiding has a universal form described 
by a single parameter, which we call the "guiding power". For short-period structures, 
d a, and for sparse structures, d <C a, it can be calculated in terms of the polarizability 
of the patterned dielectric core via a self-consistent procedure to include local field effects. 

We discuss the waveguiding in highly inhomogeneous structures, such as planar regular 
arrangements of nearly overlapping cylindrical pores or high-index cylindrical rods with 
large spacings. For thin slabs, the local field effects are different from those in an infinitely 
extended 2D PC, primarily because of the short-range dipolar interaction between the finite- 
height cylinders or spheres. Proper inclusion of the local field effects in the low-frequency 
region enables a perturbative approach with fast convergence. 

In the design of far-infrared and terahertz semiconductor lasers the optical confinement is 
an important issue. It poses a severe problem for quantum cascade lasers where intersubband 
radiative transitions require transverse magnetic (TM) polarization of emitted waves [12] . 
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With this polarization in the long wavelength region, the so-called modal confinement factor 
Ttm is known to be small for any contrast of the dielectric constants between the core and 
the claddings 3- The reason for the small values of T TM is the reduced electric field in 
the high-index core due to the boundary conditions at the core layer surface. The decrease of 
the electric field devalues traditional attempts to improve waveguiding by choosing cladding 
layers with lower index. 

In this paper, we propose to enhance the confinement by using penetration of the electric 
field into a patterned core. The model core under consideration comprises additional cylinder 
inclusions of high dielectric constant compared to that of the claddings. The filling factor 
for these inclusions is small, so as not to disturb processes in the active region (which is 
the remainder of the core outside the inclusions). We show that the better penetration of 
the field can result in a much stronger waveguiding (much reduced field spread outside the 
layer), as compared to the homogeneous layer. This case has an advantage of not being 
critically sensitive to the composition of the claddings. 

An alternative way to improve the guiding of TM waves is to use a porous PC slab as 
the active layer. We show that the confinement of TM modes in a structure with a high 
index contrast (typical for silicon-on-insulator devices) can be enlarged if the active layer 
has a patterned structure with pores. Even though the average index of the porous core is 
reduced, this is more than compensated by better penetration of the electric field into the 
structure, so that the wave confinement is ultimately enhanced. 

Patterned structures with the enhanced confinement can be advantageously used in the 
laser design to minimize the losses from free-carrier absorption and reduce the threshold 
current. 

An unusually strong guiding of TM waves in a layer of cylinders was previously observed 



in numerical studies of waveguiding by PC slab structures la, llg] but the effect had not 
been properly recognized or explained. 



II. WEAKLY GUIDED WAVES IN LATERALLY UNIFORM WAVEGUIDES 

In order to introduce the concept of weak waveguiding, let us first consider the propaga- 
tion of an electromagnetic wave along a laterally uniform dielectric waveguide with a core 
layer. Let the index profile depend only on z, approaching at large \z\ the (background) 
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dielectric constant e& of the cladding layers. 

Two wave polarizations are distinguished by the field orientation relative to the structure 
symmetry plane. Consider the case of a TE mode, when the electric field has only an in- 
plane y component and is strictly perpendicular to the wave propagation direction x. The 
wave equation for the electric field E y = E y (z) exp(iqx) is of the form 

^Ey(z) = [q 2 - k 2 e(z)]E y , (1) 

where q is the 2D wave vector in the plane of the waveguide. Let us integrate Eq. (JTJ 
between —z\ and Zi, that is over the region where the permittivity is variable. At \z\ > Z\ 
the solution of the Eq. (JTJ has the form E v {z) = E(0) exp(±«^) and we get 

-2kE v (0)= I 1 [q 2 - k 2 e(z)]E y (z)dz . (2) 

J-Zl 

In the limit of weak waveguiding, kz\ <C 1, the field E y (z) is a slowly varying function across 
the entire layer. Therefore it can be replaced in the integral by a constant value taken, e.g., 
at z = 0. Outside the guiding layer we have q 2 = Ebk\ + n 2 and hence in Eq. (j2j) we can 
take q 2 = ebk" 2 as the zeroth-order approximation. This yields 

K= \ tbkl SZ^~ i ) dz - (3) 

We have replaced the limits of integration by ±oo, since the region where e(z) = does not 
contribute to the integral in Eq. (J3]). 

Similar arguments can be used to consider the waveguiding of TM waves, when the only 
non- vanishing component of the magnetic field is H y . For this case, in the long wavelength 
limit, the quantity that remains a smoothly varying function across the layer is the normal 
component of the electric displacement vector D z . Integrating the wave equation for D z we 
obtain 

1 r 2 / I , € b 



s> kl L\ l -t)r (4) 

For both polarizations the dispersion relation for the guided wave is of the form 

e&^o = I 2 ~ ( 5 ) 

where fco = uj/c is the frequency parameter. Parameter k describing the exponential decay 
of the wave away from the core, exp(— kz), depends on frequency. Since k oc k\ the spectrum 
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of the guided wave in the long-wavelength limit has a universal character. It is convenient 
to introduce another parameter g, 

K = 2 e b k o9, ( 6 ) 

which we shall call the "guiding power" of the high-index core. The value of g defined by 
Eq. is owing to the fact that it is frequency-independent in the weak guiding limit, 
nd <C 1 [13]. So long as g is constant, Eqs. (El E]) define a universal dispersion relation for 
the guided modes of a three-layer dielectric waveguide of core thickness d. 
According to Eqs. ([3]) and (TjJ, the guiding power is given by : 

9TE = L fir " l ) dz • 9TU= L{ 1 ~^)) dz - (7) 

For the simplest case of a constant permittivity (e g ) core layer, Eqs. (JZJ) reduce to 

£te = — — —d , g T M = — — —d . (8) 

In structures with a low index contrast, (e g — e&) <C e^, the values of g for both modes 
are small and close to each other. In the opposite limit, (e 9 /e&) ^> 1, the guiding power 
for the TM mode is tb/tg times weaker than #te, which can be explained by the reduced z 
component of the electric field inside the slab. 

Confinement of guided waves is usually described by the dimensionless "confinement 
factor" T (fraction of the wave intensity that flows in the high-index core). Quite generally, 
T is proportional to the guiding power. The condition r C 1 corresponds to the weak 
guiding limit. In this limit the guiding power determines both the confinement properties 
and the dispersion of waves. 



III. WAVEGUIDING BY A PC SLAB 



Consider now the electromagnetic wave propagation along a photonic crystal slab formed 
either by a lattice of holes in a core layer or by a "mirror" structure with a set of cylinder rods 
serving as the core. In the long wavelength limit A 3> a,d, where a is the PC lattice constant 
and d is the layer thickness, the field inhomogeneity is important only at short distances 
away from the slab, since the short-range components of the fields decay exponentially over 
the distances of order a. This allows us to identify the polarized waves as TE-like and 
TM-like. For A > d the weakly bound guided waves have a smooth exponential decay of 
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the wave field away from the core, E z = E Zj out exp(— kz) with nd <C 1. The relationship 
between k and q can be obtained either by integrating the wave equation using the weak- 
guiding approximation (as in the preceding section) with averaging in the lateral plane, or 
by using the effective media approach. The latter corresponds to replacing the PC layer by 



a homogeneous slab with an effective (anisotropic) dielectric constant [T7| . Equations (0 E]) 
remain valid in both cases, but the guiding power g must now be evaluated taking account 
of the polarizability of the laterally inhomogeneous guiding layer. 

Consider a core comprising a set of dielectric cylinders of radius r with the dielectric 
constant e cy i and height d 3> r. The cylinders are spaced apart with a lattice constant a and 
the core medium between the cylinders is assumed to have the same permittivity e b as the 
background. To compare the results with those of full-scale calculations, we further consider 



a square lattice of cylinders and the structure parameters close to those studied in [15|, [16]. 

Importantly, the polarizability of a single cylinder is highly anisotropic. For a sufficiently 
elongated cylinder, d/r > 5, the polarization vector inside the cylinder is homogeneous and 
equal to that of an ellipsoid with a high axes length ratio. Besides, one must allow for depo- 
larization effects. We take them into account approximately, by replacing the cylinders by 
prolate ellipsoids with the same diameter and the same volume, so that the axes length ratio 
of the ellipsoid R e = 3/4(d/r). The dipole moment of a single cylinder in an external electric 
field equals P z = Va z E z , where V = nr 2 d is the cylinder volume and the polarizability of 
a single cylinder a z is given by 

1 e cy i — €b . . 

47r e b + (ecyi - e b )n z 



Here n z is the depolarization factor [20(. Even though in our case n z <C 1, the product 
e cy i n z cannot be neglected. The average dielectric function of the layer equals 

ei,u = e b + fe b /. cyl " " h . (10) 
e& + (e cy i - e b )n z 

The guiding power can now be calculated as in Eq. (jSJ), 

fi'TM = d(l - e b /ei t || ) , (11) 

and is given by 

9tu, cyi = d — 77- ° yl w b ^ r • (12) 

e b + [f + n z )(e cyl - e b ) 
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The spectrum of the TE mode can be calculated in a similar way, with the transverse layer 
polarizability, j_, expressed in terms of the transverse polarizability of a single cylinder, 
0L y . The result is similar to Eq. except for the depolarization factor, which should now 
be replaced by n y . For a highly elongated cylinder, n y is very close to 0.5, so that 

a y = ^ e K yl 1 7 ■ ( 13 ) 

2?r e cy i + £b 

In calculating the average dielectric function j_ of the layer one should take into account the 
local field effects (enhancement of the local field due to the field of surrounding cylinders). 
The simplest way to do this is to use the Maxwell Garnett approach, which is strictly 
applicable for dCa. In this approximation we have 

eb + e cy i + /(e cy i ~ gb) (*,s 

Cf> + Ccyl — JlCcyl ~ £b) 

The guiding power for TE waves is then given by 

STE, cyl = 2d f^~^ . (15) 

+ Ccyl — J l^cyl — CfeJ 



The calculated spectra of TM and TE modes in the "guiding power" approximation 

are shown in Fig. 1 for a square lattice structure with e cy i = 12, e g = e& = 1, d/a = 2, 

and r/a = 0.2 (n z = 0.03). The results are compared with those of full-scale numerical 

calculations [15]. Our spectra shown in Fig. 1 do not include the effect of Bragg reflection, 

which can be allowed for by using a perturbative approach If]]. With Bragg reflections 
I I 

included [181 ] . the guiding-power spectra approach the exact curves very closely. 

Note that the guiding-power spectra are described by one parameter g for each mode. We 
can determine the values of these parameters by a best fit to the exact curves. In n/a units, 
the best-fit parameters are #te = 1-2 and #tm = 3.2. These values can then be compared 
with those directly calculated from Eqs. (|12|15p . giving g TEj cy i = 1.5 and <7tm, cyl = 3.2. 
For TE waves, the Maxwell Garnett approximation apparently overestimates the local field 
effects. If we use the polarizability of separate cylinders, we find #te = 1-3, which is closer 
to the best fit. 

For the TM mode, which is of practical importance for intersubband lasers and is our main 
interest in this paper, the calculated value £/tm, cyi provides an excellent approximation. We 
regard this very good agreement as a justification for applying the guiding power approach to 
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FIG. 1: Low- frequency spectra of a photonic crystal slab composed of a square lattice of cylindrical 
rods (left panel) calculated in the "guiding power" approximation, Eqs. (0 [6]). The TE mode is 
shown by a blue solid line, the TM mode by red line. Results of full-scale numerical calculations 
are shown by open (TE) and closed (TM) dots. The dashed line indicates the light cone 
boundary. The filling factor is / = 0.125 and the cylinder's permittivity is e cy i = 12. 

calculations of the confinement factor. Note that due to the layer anisotropy, the waveguiding 
for TM waves is much stronger than that for TE waves. 

Consider now the waveguiding in structures with a porous core layer. To calculate the 
guiding power for TM waves in this case, we use a gain Eq. ( TTTT) . However, the depolarization 
effects are now very different. As is well known [20j], these effects are important when the 
high-permittivity component has a convex shape, as is the case for cylinders. In contrast, 
for the cylindrical holes, the effect of depolarization can be essentially neglected, due to the 
concave shape of the high-permittivity component. Hence ii is given by 

% || = e 3 + (e& - e g )f , (16) 

where the fraction of high-polarizability material in the core equals 1 — / for a film of porosity 
/. The resultant TM-wave guiding power is of the form 

_ H (l-/)(e g -Q.) (u] 
9 ™' P ° r " d e b + (l-f)(e 9 -e b ) ' 

To calculate the guiding power for TE waves, we again use the Maxwell Garnett approxi- 
mation, noting, however, that it is strictly applicable only for thick enough guiding layers, 
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FIG. 2: Low- frequency spectra of a photonic crystal slab formed by a dielectric layer with a 
triangular lattice of air holes with lattice constant a (left panel). Layer thickness is d = 0.6a, 
material permittivity is e„ = 12 and the hole radius is r = 0.45a. Results of the full-scale numerical 

n 

calculation [15J are shown by open (for TM wave) and closed (TE) dots; red and blue lines show 
the spectra of TM and TE waves in the "guiding power" approximation; dashed line shows the 
edges of the light cone. 

d^> r (this restriction is, of course, compatible with A ^> d required for the validity of the 
guiding power approximation). For a porous layer we find 

e b + e g + f(e b - e g ) 



e/ ' ±_ ~ 9 eb + e g -f{e b -e g ) 
whence the guiding power for TE waves is of the form 

(i- /)(<£- 3) 



(18) 



PTE, 



por 



d- 



e 9 -f(e b -e g )} ' 



(19) 



To check these results against those of full-scale numerical calculations [15J], we used 
Eqs. (0 ED to calculate the TE and TM wave spectra for a triangular lattice of pores with 
d = 0.6a, and r = 0.45a (the porosity / =0.734), see Fig. [21 The values of g in -n/a units 
that give best fit to the spectra are #te = 3.1 and #tm = 1-4. The values calculated from 
Eqs. fl!7fl9p are gxE, por = 3.4 and #tm, por = 1-4. As expected, the spectrum of the TM 
guided wave is well described with an average dielectric constant of the core layer over most 
of the Brillouin zone up to the zone boundary, where it is strongly modified by the Bragg 
reflection. The Bragg reflection can be allowed using a perturbative approach jbsl. \\\. For 
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a high-porosity layer of fairly small thickness, the spectrum of TE waves is also described 
by the guiding power approach, but the value of #te, por is overestimated in the Maxwell 
Garnett approximation. 

Note that both TM and TE waves remain weakly guided even in the extreme limit of 
sparse structures, a 3> d, so long as A > a fis| . 



IV. CONFINEMENT FACTOR OF TM WAVES BY A PC SLAB 

The modal gain in a three-layer slab waveguide with an active core layer (ACL) of thick- 
ness d can be expressed as a product of the material gain, Gacl = k n"/ €acl and the 
dimensionless "optical confinement factor", V. For weakly guided TE waves with a smooth 
variation of the electric field of the wave across the layer, the confinement factor is just 
proportional to the guiding power, Tte = nd = e\,k\dg^j2. For TM waves the confinement 
is influenced by the weakenin g of the electric field in the high-permittivity core. Generally, 
it can be written in the form [l4j 

r ™ - r Eidz ' (20) 

where the integral in the numerator is taken over the active region. For a homogeneous 
core layer one can calculate the confinement factor explicitly. In the case of weak guiding, 
nd <C 1, using Eq. ( 1201) and the boundary condition for the normal component of the field 
at the core edge planes, we find 

_ (k d) 2 ( e b \ 3 ( ^ e b 

1 TM — 7T~ 



e 9/ \ e 9 



i - - K • (2i) 



The value of Ttm is smaller than Tte by a factor (eb/e g ) 3 . For a given core composition (i.e. 
for a fixed value of e g ) and as a function of the cladding-layers index, Ttm has a maximum 
value, 

Ttm, max « 0.05 e g {k d) 2 , (22) 

achieved when the ratio of the dielectric constants is €b/e g = 3/4. Note that r T M, max is about 
ten times smaller than the confinement of the TE mode achievable in the same structure. 
For a fixed value of e g , the factor Tte does not have a maximum, except at e& = 1, where 
Tte ~ 0.5 e g (kod) 2 ~ 10 Ttm, max- The difference in confinement factors is negligible in the 
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FIG. 3: Normalized confinement factor for a PC slab with a square lattice of high-polarizability 
cylinders embedded in the core layer plotted as a function of / for several values of the ratio 
R = e cyiAg, for e g = 12 and e b = 11. 

structures with very small index contrast and thus very small confinement factors for both 



waves 



In the preceding section we showed that the TM waveguiding can be enhanced and even 
can exceed that for the TE mode by incorporating in the core layer cylindrical rods of high 
polarizability. We shall now show that the mode confinement in the active part of the core 
layer can be enhanced as well. 

Consider the waveguiding in a patterned structure with initially small difference between 
the dielectric constant of the active layer (denoted by e g ) and that of the cladding back- 



ground, €b- Smal 
active region 



values of (e g — e&) are typical for quantum cascade lasers with a multilayer 
12|. To enhance the waveguiding, we incorporate in the structure a set of 
cylindrical rods of radius r, lattice constant a, and dielectric constant e cy \. (Although we 
speak of the "lattice constant", the periodicity of rods is of no importance here and the 
result can be expressed in terms of their fill factor /.) For a waveguide with a patterned 
core layer, we first use the average dielectric constant to estimate the average field (which is 
somewhat different from the local field), and then use Eq. ( l20l) to calculate the confinement 
factor of guided TM waves 

(k d) 2 ( e b \ 2 ( \ e b 



TM 



(1 - f)e 



9 ■ 



(23) 
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FIG. 4: Normalized confinement factor as a function of e cy i for a PC slab with a square lattice 
of high-polarizability cylinders embedded in a core layer (e g = 12) for several values of the filling 
factor. The background permittivity outside the core layer is e& = 11. 

It can be seen from Eq. (1231) that in structures with a high ratio R = e cy i/ e&, the enhancement 
of waveguiding with / due to the increasing factor (1 — e^/e^ m) can overwhelm at small / 
both the decrease in (1 — /) and the decreasing factor (e^/e^ \\) 2 that describes reduction of 
the electric field in the active layer. 

The variation of the confinement factor, calculated with Eq. (1101) [modified to include a 
term proportional to (e 9 — e&)] and Eq. (123|) . is plotted in Fig. [31 The increase of confinement 
with / is mainly due to the guiding power enhanced by the better polarizability of cylinders, 
hence it is most effective at small /. It can be seen from Eqs. (|23|) and ({TO]) that the increase 
of confinement with / takes place only provided e g > 3(e g — e b ) + 1 and e cy i > 2e g . Note, 
that the increase of ii with / towards the optimum value ii opt should be evaluated at 
fixed eb and e g (we have used e b = 11 and e g = 12). This gives ^ op t = \fl$)€b- For 
e cy i 3> e 9 the optimum filling factor is small, /< 1, and the factor 1 — / is close to unity. In 
this most favorable case, the resulting confinement factor is larger than that for the initial 
homogeneous structure and approaches the optimal value given by Eq. ( |22l) . 

The dependence of the confinement factor on the cylinder dielectric constant for / = 0.05 
is shown in Fig. HI In terms of the cylinder polarizability a z (Eq. [9]), the enhancement of 
confinement is well pronounced for fa z > 20, which can be achieved in the low-wavelength 
limit by inclusions of extremely highly-polarizable materials. 
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FIG. 5: Porosity dependence of the guiding power #tm (in units of d, left panel) and the normalized 
confinement factor Ttm (right panel) for a PC slab with a porous core layer. Both parameters are 
plotted as functions of the porosity / for several values of R = e g /e^. 

Next we consider the waveguiding in structures with a porous core layer. As seen from 
Eqs. f fT7l) . the guiding power of a porous active layer is smaller than that of a homogeneous 
layer for any ratio R = e g /eb and is decreasing with /. However, for large R the decrease is 
small, while the rise of the electric field with the layer porosity provides an enhancement of 
the TM wave confinement. 

Variations of the guiding power and the confinement factor with the porosity / for several 
values of the ratio R are shown in Fig. [51 For a given R, the maximum value of Ttm is 
achieved at 

/max = 7 • (24) 



R-l 

Hence, replacing the homogeneous layer by a porous layer can increase the modal gain only 
using materials with R > 3. The maximum value of Ttm is 

4i? 3 „ 



TM, max 



horn 



(25) 



27 (R- l) 2 

where T^m is the confinement factor for a homogeneous core layer of the same R. As seen 
from Eq. fl25l) . the increase of confinement by porosity can be tangible only for structures 
with large R. For Si/SiC^ structures with R = 6 the maximal confinement enhancement is 
achieved at f — 0.6, but it is not very large (1.28); for R = 12 the enhancement is 2.1. 
We remark that although the numerical values above can be determined more accurately 
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with more elaborate numerical calculations, the increase of confinement by patterning ini- 
tially homogeneous layers is an exact result for waveguide structures with a high enough 
index contrast between the core and the claddings. The result is based solely on the power- 
law / dependencies of the guiding parameter g TM and the field ratio inside the active layer. 

V. CONCLUSIONS 

In conclusion, we considered the polarization-dependent waveguiding of light by thin 
highly inhomogeneous slabs embedded in a uniform medium. We examined exemplary slab 
structures comprising a monolayer of patterned cylindrical pores etched in a active core 
layer or a pattern of high-index dielectric rods (cylinders), embedded in the core layer. We 
demonstrated that for an optimal choice of the patterned layer structure an increase of 
optical confinement of TM wave is possible compared to a homogeneous layer. This increase 
can be achieved both by incorporating sparsely separated high-index rods which promote 
electric field penetration in the patterned structure and by using high-porosity core layers. 
Our results can be useful in the design of quantum cascade lasers. 
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